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We calculate the two- and three-loop massive operator matrix elements (OMEs) con- 
tributing to the heavy flavor Wilson coefficients of transversity. We obtain the complete 
result for the two-loop OMEs and compute the first thirteen Mellin moments at three- 
loop order. As a by-product of the calculation, the moments A*' = 1 to 13 of the 
complete two-loop and the Tf-part of the three-loop transversity anomalous dimension 
are obtained. 



1 Framework 



The transversity distribution belongs to the three twist-2 parton distribution functions 
(PDFs), together with those for unpolarized and polarized deep-inelastic scattering. It 
is a flavor non-singlet, chiral-odd distribution and can be measured in semi-inclusive deep- 
inelastic scattering (SIDIS) and via the polarized Drell-Yan process. Different experi- 
ments perform transversity measurements at the moment, cf. Refs. [2]. Recently, a first 
phenomenological parameterization has been given for the transversity up- and down-quark 
distributions in Ref . [3] , the moments of which are in qualitative agreement with first lattice 
calculations [4]. 

For semi- inclusive deeply inelastic charged lepton-nucleon scattering IN I'h + X the 
scattering cross section is given by 



SIDIS 



dxdydz 




yf] Faix,Q^)Da{z,Q'' 



-{1 - y)\S^\\Sh^\cos{^s + cl^sj ATFa{x,Q^)ATDa{z,Q^)} , (1) 



after the P?i_L-integration has been performed, [1]. We consider, for definiteness, only scatter- 
ing cross sections free of k_L- effects to refer to twist-2 quantities, x and y denote the Bjorken 
variables, z the fragmentation variable, = —q^ the space-like 4-momentum transfer, aem 
the fine structure constant, Ca the quark charge, and s the cms energy squared. S_l and 
Sh± are the transverse spin vectors of the incoming nucleon and the measured hadron h. 
Faiz^Q^), A.TFa[z,Q'^) and D a{z , Q"^) , At D a{z , Q'^) denote the unpolarized and transver- 
sity structure- and fragmentation functions, respectively. The angles 4'S,Sh ^'^s measured 
in the plane transverse to the 7*A^ axis between the x-axis and the respective vector. In 
process ([1]) the spin of the transversely polarized hadron h has to be measured. 

The transversity distribution may also be measured in the transversely polarized Drell- 
Yan processes. In Mellin space the scattering cross section is given by, [5], 



cos(2(^)ATi?(Ar, Af ^) ■ AtC^^ [N, M^) (2) 



^For a review see Ref. [1]. 
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where N denotes the Melhn variable and (j) is the azimuthal angle of one of the final state 
leptons relative to the axis defined by the transverse polarizations. 

AtH{N, Q2) = ^ [ArqiiN, Q^)ATq2{N, Q^) + Arq.iN, Q^)ATq2{N, Q^)] 

is a combination of transversity parton distributions for the incoming light (anti-)quarks, 
and AtC^^(A^, A/^) denotes the Wilson coefficient of the Drell-Yan process, with the 
invariant mass of the produced lepton pair. 

Like in the case of unpolarized and polarized deep-inelastic processes transversity receives 
heavy flavor corrections in higher orders in QCD. These are given by the corresponding heavy 
flavor Wilson coefficients. As for other non-singlet quantities [6,7], these corrections start 
at 0(aj), with as = as/(47r). In SIDIS one can tag QQ-production in the same way as in 
the deep-inelastic process, [8]. A measurement is possible in high luminosity experiments. 
In the Drell-Yan process, on the other hand, heavy flavor contributions emerge inclusively 
since there the final-state /+Z~-pairs are measured in the first place. The calculation of the 
heavy quark Wilson coefficients for ^ proceeds in the same way as in unpolarized 
and polarized deep-inelastic scattering [6,7,9,10] 

The complete Wilson coefficients for transversity can be decomposed into a light- and a 
heavy quark contribution 

As shown in [6] , the heavy quark Wilson coefficient for hard processes factorizes into the light 
quark Wilson coefficients and the massive operator matrix element at large enough 

scales ^ rn? . We apply this to the heavy flavor Wilson coefficient for transversity i?™ 
„Ta / m^X „TR.iight / ,TR ( 

= al [a.<>^-^™(7V, + 1) + ATCf\N,)] + a? [ArAf^'^'-^^iN, + 1) 

+ AT4',>™(iV, + l)®ATCW(iV, + l)+Cf' (TV,)] . (4) 

The aim of this article is to present a computation of the renormalized two- and three-loop 
heavy-flavor operator matrix elements contributing to transversity. Details of the calculation 
are given in Ref. [11]. The operator matrix element ((j'|0'^^''^^ jg) is given by a two-point 
Green's function containing a closed loop of a heavy quark Q and external massless quarks 
q. The local operator is given by 



2 



trace terms , (5) 



cf. [12]. Here S denotes symmetrization of the Lorentz indices, ti^^ = {i/2) [7^7,^ — 71^7^], 
and is the covariant derivative. The Green's function in _D = 4 -f e dimensions obeys the 
following vector decomposition 

/^y,TR,NS 7 / i^NS.TR i v r [a u« 3TR,NS at\ , AM 

G;:,,q' = JN{q\Op^^.^^^ Jq) = S,,'>^Apa''PA^^Q iy"^' ' j ' 

+ C2P^ + cs-f^'p + C4/^pA>' + C5^fp^^ (A • pf-' (6) 
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contracting the OME with a source term Jm = . . . A^" , with A^ = 0, with p the parton 
momentum. It determines the un-renormahzed massive OME 



A 



TR,NS 

'qq-Q 



AN, (A • pf+' (D - 2) 



Tr 



-A-pTr V^G^'^q' 



i A • p Tr 



" uAii.TR.NS 



(7) 



A total of 129 diagrams contribute, which were generated using QGRAF [13]. These were 

projected onto q , cf. [9], using codes written in FORM [14]. After tensor reduction, 
the loop integrals are of the tadpole-type, since the single external quark is on-shell and 
masslcss. The integrals were then evaluated using MATAD [15]. The renormalization of the 
OMEs is described in Ref. [9]. 

After mass- and charge renormalization one obtains the massive OMEs in the on-mass- 
shell scheme, cf. [9], 



. 4(2),NS,MS 



. -(3),NS,MS 



a (0),TR 
pQ,Qrqq' ^^2 



- (1),TR 



In 



(0),TR 
(2),TR PO.Qlqq 



(8) 



^(/3o + 2/3o,g)ln^(^)+i(27«^™/3o. 



-27W'™(^o + /3o,q)+/?i,q7^"/ 




^(2),TR 
Iqq 



-( 



Aa 



(2),TR 

qq,Q 



(0),TR r\ > yLj^^[\.r} > 

_(2),TR.^ „ X ^,(0)o(2) Iqq PoPO.Q(,3 Iqq PG,Q(,2 

qq^Q (Po + Po,q)-7VPi,q Q 1 



„(i).TR, 



-U9Arr7(l)« ^(0),TR , r^(0)r,(l),TR , r,r (-1) (2),TR (3),TR 



+4:a 



at 2- and 3-loops. Here, Cfe denotes the Riemann ^-function at "fqq''^^ are the transversity 
anomalous dimensions for / = 0, 1, 2 in LO [16], NLO [5, 17], and NNLO [18], with f{Nf) = 
f{Nf + 1) — f{Nf). For the other quantities we refer to [9]. The new terms being calculated 



are a 



(2),TR 

qq,' 



(3),TR 



(2),TR, 



2 Results 

2.1 Massive Operator Matrix Elements 



At 0{aj.) the massive operator matrix elements for transversity AtA^^^q^'^^ are obtained 
for general values of N, cf. Eq. ([8|). The im-rcnormalized OME is computed to 0{e) to 

also extract the functions a™n'^\N). The new terms at 2-loops are aJfVi^'' and aJfVi^\ cf- 

qq^w qq^w qq^w 
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Eqs. 



a 



TR,(2) 



(N) = CfTf 



-S3 



40 



S2 



"224 
'27 



+ T7C2 



Si + 2C2 



(24 + 73 TV + 73 
18iV(l + iV) 



(10) 



656 20 8 ■ 



5i 



112 4 
^+3^' 



4 12 

^4 + 77 C2 + TT Cs + 

3 6 3 



-144-48iV + 757iV2 



S2- 
1034 iV3 



216iV2 (1 + iV)^ 



20 

-517N-^ 



(11) 



with Sk = Sk{N) the single harmonic sums. 

At 0{al) the moments iV = 1 to 13 were computed for the massive OMEs, as e.g. 



At ^-^(13) 



'1751446 



110565 
20032048197492631 

-C-F - 



7005784 , A , 3/m2 

^--T2162T5^-(^^ + ^0l-(^ 



2193567563187000 8027019000 
/ 1705832327329042449983 7005784 



137401473299^ 93611152819 



A 263491335690022440000 45045 
7005784 ^ \ „ 458114791076413771 



3652293645 'J V /i^ 
/ 3385454488248191237 
V 65807026895610000 



-( 



45045 
7005784 



C3)Ca 



6580702689561000 



3648645 



3502892 ^ 
Bi + C4 



81735983092 



C3 



135135 15015 243486243 

55376278299522733837425052493\ / 3502892 

122080805651901196900800000 / ^ ^ V 135135 
4061479439 3486896974743882556775647\ 

12162150 ^ 12935029206601101600000 / ^ 
/ 279922752632160355860697 56046272 \ 

V 3557133031815302940000 ^ 1216215 / ^ ^ 

/ 291526550302760070155303 14011568 



3502892 
15015 ^ 



V 7114266063630605880000 



173745 



(12) 



where 



-4C2ln2 (2) + -ln4(2) 



13 



C4 + 16Li4 



Like for the massive OMEs in case of unpolarized deep-inelastic scattering, the structure 



of At a 



(3),NS,MS 



[N) is widely known for general values of N , except the contributions due 



to the finite part o/^qQ and the 3-loop anomalous dimension Yqq {^)- One notices the 
cancellation of all C2 terms in A!^qq Q^'^^{N) after renormalization. 
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2.2 Anomalous Dimensions 



The transversity anomalous dimension is given by 

oo 

■(iV,a,) = ^<7^;)'™(iV). (13) 



yTR, 
Iqq 



From Eq. ^ one may determine the complete 2-loop anomalous dimension [5,17] and the 
Tp-part of the 3-loop anomalous dimension [18]. We agree with the results given in [5,17] and 
eonfirm the Tp-contributions for the moments A'^ = 1 to 8 given in Refs. [18]. Furthermore, 
we obtain 'jqq^'™' = 7gq^'™(-^/ + 1) ^ lqV''^^i^f) newly for = 9 to 13, as e.g. 



^(3),TR ^ 13) ^ _ CfTf 



36713319015407141570017 14011568 
131745667845011220000 ^ 45045 ^ ^ ~ '^''"^ 



66409807459266571^ 6571493644375020121 
3290351344780500 ^ ^'^ 65807026895610000 



2.3 A Remark on the Soffer Bound 



If the Soffer inequality [19] 

\ATf{x, Q')\ < \ [f{x, + A/(2;, Q2)] (14) 

holds for the non-perturbative PDFs in Eq. (|14p one may check its generalization from 
fi Fi for the corresponding structure functions. This includes the non-singlet evolution 
operator (Eq. (6), Ref. [20]) and the heavy flavor Wilson coefficient. At perturbative scales, 
it holds for the evolution operator [11], generalizing a result from [5] for the moments = 1 
to 13 at 3-loops. For the heavy quark Wilson coefficients in SIDIS we only know the massive 
OMEs so far. As shown in Ref. [11], a final conclusion can only be drawn knowing the yet 
undetermined massless Wilson coefficients. Here the difference [A^q^Q ~ ^™q] {^) of t^ie 
massive OMEs, shows a sign change to negative values for /m^ in the physical range. For 
large scales Q^/rri^ ^ 1 positive values are obtained. 
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